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Abstract 

Finite-size scaling (FSS) of the five-dimensional {d = 5) Ising model is investigated numerically. 
Because of the hyperscaling violation in d > 4, FSS of the d = 5 Ising model no longer obeys 
the conventional scaling relation. Rather, it is expected that the FSS behavior depends on the 
geometry of the embedding space (boundary condition). In this paper, we consider the cylindrical 
geometry, and explore its influence on the correlation length £ = L n f{eL y t , HL y h) with system 
size L, reduced temperature e, and magnetic field H; the indices, yf h , and ft, characterize FSS. For 
that purpose, we employed the transfer-matrix method with Novotny's technique, which enables 
us to treat an arbitrary (integral) number of spins N = 8, 10, . . . , 28; note that conventionally, 
N is restricted in N(= L d_1 ) = 16,81,256,... . As a result, we estimate the scaling indices as 

= 1.40(15), y\ = 2.8(2), and y* h = 4.3(1). Additionally, under postulating (7 = 4/3, we arrive at 
y\ = 2.67(10) and y^ = 4.0(2). These indices differ from the naively expected ones, Q = 1, yj! = 2 
and y*^ = 3. Rather, our data support the generic formulas, f2 = (d — l)/3, y* t = 2(d — l)/3, and 
yt = d — 1, advocated for the cylindrical geometry in d > 4. 
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I. INTRODUCTION 



The criticality of the Ising model above the upper critical dimension (d > 4) belongs to the 
mean-field universality class. However, the finite-size effect, namely, the finite-size-scaling 
behavior, is not quite universal because of the violation of the hyperscaling in d > 4. (From a 
renormalization-group viewpoint, this peculiarity is attributed to the presence of "dangerous 
irrelevant variable" [l[ [2], y|.) Actually, it is expected that the embedding geometry of 
;he system (boundary condition) would affect the finite-size-scaling behavior for d > 4 



Recently, Jones and Young performed an extensive Monte Carlo simulation for the d = 5 
Ising model embedded in the (periodic) hypercubic geometry fl^ . They calculated the 

131 ] . (Because the calculation of £ requires a 



correlation length £ with Kim's technique 



computational effort, Binder's cumulant ratio rather than £ has been studied extensively so 
far [3, l8|, [9| . ) They found that the correlation length £ obeys the scaling relation, 

£ = L Q /(eL^), (1) 

with the linear dimension L and the reduced temperature e. They found that the scaling 

n nti 

indices are in good agreement with the theoretical prediction |4J, [5|, 16[ , 



Q = 5/4 and y* = 5/2. 



(2) 



Notably enough, their indices exclude the naively expected values, Q = 1 and yl = 2. 
That is, the correlation length £ ~ L Q exceeds the system size L at the critical point e = 
as L — > 00. (In other words, the spin-wave excitation costs very little energy for large L.) 
Such a peculiarity should be attributed to the violation of the conventional scaling relation 
(hyperscaling) in d > 4. It would be intriguing that the above formula is cast into the 



expression 



7], 



Z = lf(d 2 ), (3) 

with the replacements I = L n and yl = 2Q. The expression is now reminiscent of the 
conventional scaling relation expected for the mean-field universality class. Namely, the vio- 
lation of hyperscaling is reconciled (absorbed) by the replacements, and the scaling indices, 
Q and y^, characterize the anomaly quantitatively. So far, numerous considerations have 
been made 
very efficiently 



12J for the hypercubic geometry, where the Monte Carlo method works 



In this paper, we investigate the d = 5 Ising model embedded in the cylindrical geometry; 
namely, we consider a system with infinite system size along a particular direction. Clearly, 
the transfer-matrix method well suits exploiting such a geometry. However, practically, the 
transfer-matrix method does not apply very well in large dimensions d > 3 because of its 
severe limitation as to the available system sizes. 
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In order to resolve this limitation, we implemented Novotny's technique Q, [l5|, [16, 



n 



181 ]. which enables us to treat an arbitrary number of system sizes N = 8, 10, ... , 28; 



here, the system size N denotes the number of constituent spins within a unit of the 
transfer matrix (Fig. [T]). Note that conventionally, the system size is restricted in 
N(= = 16, 81, 256, . . . , which soon exceeds the limit of available computer resources. 

Such an arbitrariness allows us to treat a variety of system sizes, and manage system- 
atic finite-size-scaling analysis. Actually, with the scaling analysis, we obtained the indices 
= 1.40(15), y* t = 2.8(2), and y* h = 4.3(1). Moreover, postulating ft = 4/3, we obtained 
yl = 2.67(10) and yl = 4.0(2). (Here, the exponent y\ denotes the scaling dimension of the 
magnetic field. Our scaling relation, Eq. (Tl3"l) . incorporates the magnetic field H and the 
corresponding scaling index y* h .) Obviously, our results exclude the naively expected ones, 
Vt = 1, yl = 2, and yl = 3. Rather, our data seem to support the generic formulas [4, 0, 

Q = ^y^, y* t = 2-^, and y* h = d-l, (4) 

advocated for the cylindrical geometry in d > 4. Actually, our data deviate from the above 
mentioned values for the hypercubic geometry, Eq. ([2]), indicating that the embedding 
geometry is indeed influential upon the finite-size scaling. 

In fairness, its has to be mentioned that Novotny obtained Q + y[ w d — 1 Q, [l3] ■ He 
postulated Q = (d — l)/3 in order to fix the location of the critical point. In this paper, we do 
not rely on any propositions, and estimate the indices independently. For that purpose, we 
calculated the cumulant ratio to get information on the critical point. Moreover, we treated 



the system sizes up to N = 28, which is substantially larger than that of Ref. [17j, N < 13. 
Here, we made use of an equivalence between the d = 5 Ising model and the quantum d = 4 
Ising model; the latter is computationally less demanding. We also eliminated insystematic 
finite-size corrections by tuning extended coupling constants; see our Hamiltonian (jHJ). In 
this respect, the motivation of the present research is well directed to methodology. 

The rest of this paper is organized as follows. In Sec. [Til we explain our simulation 
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scheme in detail. In Sec. 1111} we manage the finite-size-scaling analyses of the simulation 
data. In the last section, we present the summary and discussions. 

II. NUMERICAL METHOD 

In this section, we explain the numerical method. First, we argue the reduction of the 
d — 5 Ising model to the d = 4 quantum transverse-field Ising model. The reduced (quantum 
mechanical) model is much easier to treat numerically. Second, we explicate Novotny's 
transfer-matrix method. We place an emphasis how we extended his formalism to adopt the 
quantum-mechanical interaction. 

A. Reduction of the classical d = 5 Ising model to the d = 4 quantum counterpart 

The c?-dimensional Ising model reduces to the (d — 1) -dimensional transverse-field Ising 
model; in general, the d-dimens.ona. classical 8y8t e,n has its {d - ^-dimensional quantum 
counterpart [19(. Such a reduction is based on the observation that the transfer-matrix 
direction and the (quantum) imaginary-time evolution have a close relationship. Actually, 
the quantum Hamiltonian is an infinitesimal generator of the transfer matrix. Because the 
quantum Hamiltonian contains few non-zero elements, its diagonalization requires reduced 
computational effort. The significant point is that the universality class (criticality) is main- 
tained through the mapping. 

To be specific, we consider the following d = 4 transverse-field Ising model with the 
extended interactions. The Hamiltonian is given by, 

(ij) ««» M IM] i » 

Here, the operators {cr"} denote the Pauli matrices placed at the d = 4 hypercubic lattice 
points i. The parameters V and H stand for the transverse and longitudinal magnetic fields, 

respectively. The summations, Yutij)i S((ij))> Stop anc ^ S [[#]]' run over an possible nearest- 
neighbor pairs, the next-nearest-neighbor (plaquette diagonal) pairs, the third-neighbor 
pairs, and the fourth- neighbor pairs, respectively. The parameters { Ji} {i = 1, 2, 3, 4) are the 
corresponding coupling constants. Hereafter, we regard J\ as a unit of energy (Ji = 1), and 
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tune the remaining coupling constants ^2,3,4 so as to eliminate the insystematic finite-size 
errors; see Sec. Mil 

We simulate the above d = 4 quantum Ising model with the numerical diagonalization 
method. The diagonalization of such a high-dimensional system requires huge computer- 
memory space. In fact, the number of spins constituting the d = 4 cluster increases very 
rapidly as N = 16, 81, 256, . . . , overwhelming the available computer resources. In the next 
section, we resolve this difficulty through resorting to Novotny's transfer-matrix formalism. 



B. Constructions of the Hamiltonian matrix elements 



In this section, we present an explicit representation for the Hamiltonian ((Sj). We make use 



of Novotny's method 14j, which enables us to treat an arbitrary number of spins constituting 
a unit of the transfer matrix. Novotny formulated the idea for the classical Ising model. In 
this paper, we show that his idea is applicable to the quantum Ising model, Eq. (jSJ), as well. 

Before we commence a detailed discussion, we explain the basic idea of Novotny's method. 
In Fig. [TJ we presented a schematic drawing of a unit of the transfer matrix for the Ising 
model in d — 3 (rather than d — 5 for the sake of simplicity). Because the cross-section 
of the d = 3-dimensional bar is d = 2-dimensional, the transfer-matrix unit should have a 
d = 2-dimensional structure. However, in Fig. [fl the spins {<7j} (o"j = ±1, i = 1, 2, . . . , N) 
constitute ad— 1-dimensional (zig-zag) structure. This feature is essential for us to con- 
struct the transfer-matrix unit with an arbitrary (integral) number of spins. The dimen- 
sionality is lifted to d — 2 effectively by the long-range interactions over the yN th-neighbor 
pairs; owing to the long-range interaction, the N spins constitute a \/N x \/N rectangu- 
lar network. (The significant point is that the number y/~N is not necessarily an integral 
nor rational number.) Similarly, the bridge over (s/N ± l)th neighbor pairs introduces the 
next-nearest-neighbor (plaquette diagonal) couping with respect to the y/N x \/jV network. 

We apply this idea to the case of the d = 4 quantum system. To begin with, we set 
up the Hilbert-space bases {|eri, ct 2 , . . . , &n)} (°"i — ±1) f° r the quantum spins {erf} (i = 
1, 2, . . . , N). These bases diagonalize the operator <rj; namely, 

^lfo}> = *ilM>, ( 6 ) 

holds. 
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We consider the one- and two-body interactions separately. Namely, we decompose the 
Hamiltonian (jSJ) into two sectors, 

n = n {2 \{j l }) + n (1 \T,H). (7) 

The component originates from the spin-spin interaction, which depends on the ex- 
change couplings {Ji}- On the other hand, the contribution H.^ comes from the single-spin 
terms, depending on the magnetic fields, T and H. 

First, we consider TiS 2 ' . This component concerns the mutual connectivity among the 
N spins, and we apply Novotny's idea to represent the matrix elements. We propose the 
following expression, 

^ (2) = y E H ^ + Y E E 

a=±«,ogA a,/3(zA a=±a B=±3 

E EE X>(a + ^ + 7) 

a,/3,7£A o=±a p = ±p 7=±7 

+ yE E E E H(a + P + l + h (8) 

a=±l p = ± N l/4 ^ = ± N l/2 g =±N 3/4 

Here, the set A consists of the elements, A = {1, iV 1 / 4 , iV 1//2 , iV 3 / 4 }. The component H(y) 
denotes the vth-neighbor interaction for the iV-spin alignment, 

Hw.UnM = {{*i}\H(y)\{n}) = {M\TP°\{n}), (9) 
with the exchange-interaction matrix, 

N 

({<T t }\T\{Ti}) = -J2^r k , (10) 

k=l 

and the translational operator P satisfying, 

P\Wi}) = I ton}), (ii) 

under the periodic boundary condition. The insertion of P v beside the T operation is a key 
element to introduce the coupling over the v th-neighbor pairs. The denominator 2 in Eq. 
dHJ) compensates the duplicated sum. 

Let us explain the meaning of the above formula, Eq. (181) . more in detail. As shown 
in Fig. [H in the case of d — 2, we made bridges over iV^th- neighbor pairs to lift up the 
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dimensionality to d = 2 effectively. In the case of d = 4, by analogy, we introduce the 
interaction distances such as v — l,^ 1 / 4 ,^ 1 / 2 and iV 3 / 4 . The first term in Eq. flS} thus 
represents the nearest-neighbor interactions (with respect to the d = 4-dimensional cluster). 
Similarly, the remaining terms introduce the long-range interactions. For example, the com- 
ponent H{1 + iV 1//4 ) introduces the next-nearest-neighbor (plaquette diagonal) interaction. 
We emphasize that the idea of Novotny is readily applicable to the quantum model. [In 
short, our (quantum mechanical) formulation is additive. On the contrary, Novotny's orig- 
inal formulation is multiplicative, because his original formulation concerns the Boltzmann 
weight rather than the Hamiltonian itself.] 

Lastly, we consider the one-body part . The matrix element is given by the formula, 

KL{n} = <{^}l^ (1) l^}>- (12) 

The expression is quite standard, because the component H.^ simply concerns the individual 
spins, and has nothing to do with the connectivity among them. 

The above formulas complete our basis to simulate the quantum Hamiltonian ([5]) numer- 
ically. In the next section, we perform the numerical simulation for AT = 8, 10, ... , 28. 

III. NUMERICAL RESULTS 

In Sec. HIl we set up an explicit expression for the Hamiltonian (jSJ); see Eqs. (jSJ) and 
(fl2"|) . In this section, we diagonalize the Hamiltonian for iV = 8, 10, . . . , 28 with the Lanczos 
algorithm. We calculated the first-excitation energy gap AE (rather than £). The scaling 
relation for AE is given by, 

AE = L~ n g(eL y * , HL y * h ), (13) 

because AE ~ l/£ holds. (As compared to Eq. p]), our scaling relation is extended to 
include the magnetic field H as well as the corresponding scaling index y^.) The reduced 
temperature e is given by e = T — T c with the critical point T c . Note that the linear dimension 
L satisfies L = N 1 ^, because the A^ spins constitute the d = 4-dimensional cluster. 
We fix the interaction parameters to, 

( Jx, J 2 , J 3 , J 4 ) = (1, 0.15, 0.05, 0.05), (14) 
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and scan the transverse magnetic field T. (We will also provide data for (Ji, J 2 , J3, J4) = 
(1,0,0,0) and (1,0.1,0.1,0.05) as a reference.) The interaction parameters, Eq. (P3J, are 
optimal in the sense that the insystematic finite-size errors are suppressed satisfactorily. 



A. Scaling behavior of Binder's cumulant ratio and the transition point 

Because the scaling relation, Eq. ([TBI , contains a number of free parameters, it is am- 
biguous to determine these parameters simultaneously. Actually, in Ref . [14], [17j , the author 
fixed Q = 4/3, to determine the index y* t . 

In this paper, we estimate the scaling indices independently without resorting to any 
postulations. For that purpose, we calculated an additional quantity, namely, Binder's 



cumulant ratio 



20|, 



to determine the location of T c . Here, the brackets (...) denote the expectation value at the 
ground state. The magnetic moment M is given by M = YliLi (J i- Because the cumulant 
ratio is dimensionless (Q = 0), it obeys a simplified scaling relation; 

U = U (eL y t , H L y * h ) . (16) 

Hence, the intersection point of the cumulant-ratio curves indicates a location of the critical 
point. The scaling relation for the cumulant ratio, Eq. ffTBT) . has been studied extensively 

nhfln 

for the d = 5 hypercubic geometry with the Monte Carlo method [7|, [8|, [9|, [12J. 



In Fig. [21 we plotted the cumulant ratio for various T and N = 8, 10, . . . , 28 with H = 
fixed; as mentioned above, we fixed the exchange-coupling constants { J,} to Eq. (1T4"|) . From 
the scale-invariant (intersection) point of the curves in Fig. [2J we observe a clear indication 
of criticality at r c w 12.5. In the subsequent analysis of Sec. IIIIB| we make use of this 
information to determine the scaling indices. 

This is a good position to address why we fixed the exchange couplings to Eq. f[T4j) . As 
a comparison, we presented the cumulant ratio for various T and N = 8, 10, . . . , 28 in Fig. 
[21 where we tentatively turn off the extended couplings «/2,3,4 = 0. Clearly, the data are 
scattered as compared to those of Fig. [2J Such data scatter obscures the onset of the phase- 
transition point, and prohibits detailed data analysis of criticality. In order to improve the 
finite-size behavior, we surveyed the parameter space {Ji}, and found that the choice ffl^l) is 
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an optimal one. Such elimination of finite-size errors has been utilized successfully in recent 
numerical studies [2ll. \22\. 



B. Critical exponent 

Provided by the information on T c (Fig. [2]), we are able to determine the scaling indices 
from the scaling relation (TT3T) . In this section, we consider the index Q. 
In Fig. HI we plotted the approximate index, 

n(Ll ' h2) ~ HMIL2) ' ( } 

for [2/(I,i + L 2 )] 3 with 8 < Nt < N 2 < 28; note that L 1>2 = N^ 2 4 holds. The parameters 
are the same as those of Fig. [2j The approximate transition point r c (Li, L 2 ) is given by the 
intersection point of the cumulant ratio for a pair of (Ni,N 2 ); namely, it satisfies 

U(N X , r c (L a , L 2 )) = U(N 2 , r c (L 1; L 2 )). (18) 

The least-squares fit to these data yields Q = 1.403(46) in the thermodynamic limit. We car- 
ried out similar data analysis for (J 2 , J 3 , J 4 ) = (0.1,0.1,0.05), and obtained Q = 1.494(21). 
As an error indicator, we accept the difference between them. As a consequence we estimate 
the index as 

Q= 1.40(15). (19) 

Let us mention a number of remarks. First, our result excludes the naively expected 
one fl — 1. Actually, the result O > 1 indicates that the correlation length L Q develops 
more rapidly than the system size L enlarges. This feature may reflect the fact that the 
spin waves cost very little energy large L. Second, our result supports the generic formula 
Q = 4/3, Eq. (j3J), advocated for the cylindrical geometry in d > 4 

Baa 

. On the contrary, 

it deviates from that of the hypercubic geometry (|2l) ; we confirm this observation in the 
following sections. Lastly, the validity of the abscissa scale (extrapolation scheme), 1/L 3 , in 
Fig. H] is not clear. In Sec IIII D\ we inquire into the validity of the extrapolation scheme. 
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C. Critical exponents yl and y* h 



and, 



In Figs. [5] and El we plotted the approximate indices, 

- m l.) + W,^ = wo/ww, (21) 

respectively, for [2/(Li + L 2 )] 3 with 8 < N\ < N 2 < 28; the parameters are the same 
as those of Fig. [2j The least-squares fit to these data yields the estimates, — Q + y\ = 
1.396(21) and — Q + 2yl = 7.198(39), in the thermodynamic limit. Similarly for (J 2 , J3, J4) = 
(0.1, 0.1, 0.05), we obtained -Q + y* t = 1.359(12) and -SI + 2y* h = 7.211(24). Consequently, 
we estimate the scaling indices as — Q + yl = 1.4(1) and — Q + 1y* h = 7.2(1). Combining 
them with Q = 1.40(15), Eq. (fl9|) . we arrive at, 

V*t = 2-8(2), (22) 

and, 

V* h = 4-3(1). (23) 

Again, our data exclude the naively expected values, yl = 2 and yl = 3. Rather, 
our estimates are comparable with the generic formula, Eq. (j4j); actually, the estimate 
yl = 2.8(2) is quite consistent with the prediction yl = 8/3, Eq. (Tjfl), whereas the result 
yl = 4.3(1) and the formula yl = 4, Eq. (TjO), are rather out of the error margin. We 
attain more satisfactory agreement between the numerical result and the formula by the 
data analysis under the assumption Q = 4/3 in the next section. On the contrary, our data 
conflict with the values, Eq. (T5]), anticipated for the hypercubic geometry. Hence, the data 
suggest that the embedding geometry is influential on the finite-size scaling above the upper 
critical dimension. We confirm this issue more in detail in the next section. 

D. Scaling indices yl and yl under the assumption $7 = 4/3 

In Sec. IIIIB| we obtained an estimate Q = 1.40(15) being in good agreement with the 



formula Q = 4/3, Eq. (jlj). In this section, we assume Q = 4/3 
remaining indices yl and yl under this hypothesis. 
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and estimate the 



In Fig. [7J we plotted the approximate index, 

Vt (Ll ' U) ~ HL1/L2) ' (24) 

for the abscissa scale [2/(Li + L 2 )] 3 with 8 < N\ < iV 2 < 28; the parameters are the same 
as those of Figs. HE Because we assumed f2 = 4/3 Q|, we are able to determine the 
approximate critical point r c (.Li, L2) from the fixed point of L 4//3 AL(L); namely, 

Li /3 AE(N X , f c (Li, La)) = L^ /3 AL(iV 2 , f c (L 1; L 2 )). (25) 

We notice that the data exhibit improved convergence to the thermodynamic limit. The 
least-squares fit to these data yields y\ = 2.671(49) in the thermodynamic limit. Similarly, 
we obtained y\ = 2.697(42) for (J 2 , J3, J4) = (0.1,0.1,0.05). Consequently, we estimate, 

y* t = 2.67(10). (26) 

This result is consistent with the above estimate y% = 2.8(2), Eq. (122]) . confirming the 
reliability of our analyses in Figs. HE It is also in good agreement with the prediction 
y* t = 8/3, Eq. (jlj). On the contrary, our estimate excludes the exponent y\ = 5/2, Eq. (j2J), 
advocated for the hypercubic geometry. 

Similarly, in Fig. El we plotted the approximate index, 

., r r , Hn(^^(iV 1 )/^^(iV 2 ))| w=fe(Li;L2) 
V ^ U) = 2 HL^) ' (2?) 

for the abscissa scale [2/(Li + L 2 )] 3 with 8 < iVi < AT 2 < 28; the parameters are the same 
as those of Figs. SB The data exhibit an appreciable systematic finite-size deviation. The 
least-squares fit to these data yields y* h = 4.021(60). Similarly, we obtained y* h = 4.148(36) 
for (J 2 , J3, J4) = (0.1,0.1,0.05). Consequently, we estimate, 

y* h = 4.0(2). (28) 

Again, the result is quite consistent with the prediction yl = 4, Eq. (j3J. In other worlds, 
this agreement indicates that the extrapolation scheme (abscissa scale) 1/L 3 is sensible. 

Let us mention a few remarks. First, the data in Figs. [7] and M exhibit suppressed 
finite-size corrections owing to the assumption Q — 4/3. This feature was observed in Refs. 
30]; where the author estimated yl reliably under Vt = (d — l)/3; see the Introduction. 
Our analysis shows that the assumption yields a reliable estimate for y^ as well as y\. Second, 

11 



our data confirm the self-consistency of our analyses performed in Figs. HHHJ Particularly, 
the data justify the extrapolation scheme with the abscissa scale 1/L 3 . As a matter of fact, 



in Ref. 



231 ]. the authors observed notable finite-size corrections to the cumulant ratio U 



obeying L d 2y h. In our case (d — 5 cylindrical geometry), the power should read d—1y* h = —3. 
Hence, our numerical data support their claim. Lastly, as to the convergence of U to the 



3, 



25 



20,123, 



28 



29 



30 



31 



32J ; it has been 



thermodynamic limit, there arose controversies 
reported that there appear unclarified finite-size corrections to U, which prohibit us to take 
reliable extrapolation to the thermodynamic limit. In this paper, we avoided the subtlety by 
eliminating finite-size errors with the (finitely-tuned) extended interactions; see Figs. [2] and 
El We consider that such a technique would be significant for the study of high-dimensional 
systems, where the available system size is restricted. 

IV. SUMMARY AND DISCUSSIONS 



We studied the finite-size-scaling behavior of the d — 5 Ising model embedded in the 
cylindrical geometry. Our aim is to see an influence of the embedding geometry (boundary 
condition) on the scaling relation, Eq. ffTB"]) ; the embedding geometry should alter the scaling 
indices £7 and y* th above d > 4 0, 5), jfj. For that purpose, we employed the transfer-matrix 
method (Sec. [II]), and implemented Novotny's technique 14| to treat a variety of system sizes 
N — 8, 10, . . . , 28. Moreover, we made use of an equivalence between the d = 5 (classical) 
Ising model and its d — 4 quantum counterpart; the latter version is computationally less 
demanding with the universality class retained. 

We analyzed the simulation data with the finite-size scaling relation, Eq. (|T3j) . and 
obtained the scaling indices as Q = 1.40(15), y\ = 2.8(2), and y* h = 4.3(1). Additionally, 
under fl = 4/3, we estimate y\ = 2.67(10) and y* h = 4.0(2). The indices exclude the naively 
expected ones, ft — 1, y* — 2, and yl = 3, reflecting the violation of hyperscaling in large 
dimensions. Clearly, our data support the generic formulas, Eq. (j3]), advocated for the 
cylindrical geometry in d > 4 0, 0, Q] . On the contrary, our data conflict with the values for 
the hypercubic geometry, Eq. ([2]). Our result demonstrates that the embedding geometry 
is indeed influential on the scaling indices. 

Lastly, let us mention a few remarks. First, we stress that the violation of hyperscaling 
above the upper critical dimension d is not necessarily an issue of pure academic interest. For 
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example, a class of long-range interaction 23j suppresses the upper critical dimension to an 



experimentally accessible regime d < 3. Second, there arose controversies 
29. 



24 



25 



26 



30 



31 



21, 



32| concerning the subdominant finite-size effect (corrections to scaling) above 



d. More specifically, the Binder-cumulant data exhibit unexpectedly slow convergence to 
the thermodynamic limit. In this paper, we avoided this subtlety by extending (tuning) the 
exchange- coupling constants to Eq. (JUJ), where we observed eliminated finite-size errors. 
Actually, from Figs. [2] and [31 we notice that the elimination was successful. Our data 
indicate that the (dominant) finite-size errors obey the power law L~ 3 as claimed in Ref. 
23)]. We consider that the elimination of finite-size errors is significant for the study of 
high-dimensional systems, where the available system size is restricted severely. 
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FIG. 1: Construction of the spin cluster for the quantum transverse-field Ising model, Eq. ([5]). 
For simplicity, we consider the case of d = 2. As indicated above, the spins constitute a d = 1- 
dimensional alignment {cj} (i = 1,2,..., iV), and the dimensionality is lifted to d = 2 by intro- 
ducing the bridges (long-range interactions) over the N^Hh and 1 + N l / 2 th neighbor pairs; these 
interactions correspond to the nearest-neighbor and the next-nearest-neighbor interactions, respec- 
tively, with respect to the d = 2 cluster. In the case of d = 4, we consider the N^Hh, N^Hh and 
iV 3 / 4 th neighbor interactions; see Eq. ([5D for detail. This idea, namely, Novotny's method, was 
first developed for the classical Ising model [l^j]. We apply this method to the quantum system. 
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FIG. 2: Binder's cumulant ratio U (|15j) is plotted for the transverse magnetic field V and the 
system sizes N = 8, 10, . . . , 28 with the fixed exchange couplings, Eq. (|14p . We observe a clear 
indication of criticality at r c ~ 12.5. Apparently, the finite-size-scaling behavior is improved as 
compared to that of Fig. [3l where we turned off the extended interactions «/2,3,4 = 0. 
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FIG. 3: Tentatively, we turned off the extended interactions (^2,3,4 = 0), and calculated the 
cumulant ratio U ([15]) for various T and N = 8, 10, . . . , 28. We notice that the data are scattered 
as compared to those in Fig. 
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FIG. 4: The approximate critical index fi(Li, L 2 ) (HZD is plotted for [2/(Li +L 2 )] 3 with 8 < N% < 
iV 2 < 28 (X*i,2 = ^i^ 4 ); the parameters are the same as those of Fig. [2l The least-squares fit to 
these data yields Vt = 1.403(46) in the thermodynamic limit L — ► 00. 
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FIG. 5: The approximate critical index L 2 ) + yt(Li, L 2 ) ((20]) is plotted for [2/(Li + L 2 )] 3 

with 8 < N\ < N2 < 28 (£1,2 = -^1 2 4 )j the parameters are the same as those of Fig. [2j The 
least-squares fit to these data yields — O + y\ = 1.396(21) in the thermodynamic limit L — > 00. 
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FIG. 6: The approximate critical index -0(L X , L 2 ) + 2y£(Li, L 2 ) <^B) is plotted for [2/(L 1 + L 2 )} 3 
with 8 < A^i < N 2 < 28 (£1,2 = ^i^ 2 )] the parameters are the same as those of Fig. [2j The 



least-squares fit to these data yields — + 2y* h = 7.198(39) in the thermodynamic limit L 



oo. 
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FIG. 7: The approximate critical index yf(L 1 ,L 2 ) (JUD is plotted for [2/(Li +L 2 )f with 8 < N\ < 
N 2 < 28 (^i ; 2 = ^i^ 4 )' the parameters are the same as those of Fig. (2) The least-squares fit to 
these data yields y\ = 2.671(49) in the thermodynamic limit L — > oo. The result is consistent with 
Eq. (|22p . confirming the reliability of our analysis. 
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FIG. 8: The approximate critical index y* h (L l7 L 2 ) is plotted for [2/(L 1 +L 2 )f with 8 < N% < 
N% < 28 (Xi,2 = )i the parameters are the same as those of Fig. [2l The least-squares fit to 
these data yields y* h = 4.021(60) in the thermodynamic limit L — ► 00. The result is consistent with 
the prediction ([!]). 
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